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The angular deficit factor in the Levi-Civita vacuum metric has been parametrized 
using a Riemann-Liouville fractional integral. This introduces a new parameter 
into the general relativistic cylinder description, the fractional index a. When the 
£>V fractional index is continued into the negative a region, new behavior is found in the 

Gott-Hiscock cylinder and in an Israel shell. 



I. INTRODUCTION 

The 1917 Levi-Civita x ' solution is a standard vacuum exterior for cylindrical matter 
distributions. The metric can be written in terms of two parameters: mass density a and 
angular deficit b 
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The interior matches add a mass/length parameter, A, to the description of the complete 
solution. Currently 3-dim cylindrical matter distributions with positive density are believed 
to provide a good description for < a < 1/2.I3~UJ The literature has tended to focus on 
the (cr, A) ranges but one additional feature in many of the interior solutions is the boundary 
location and size that varies as a function of the density. The size of the circumference 



remains finite even as the density increases and a coordinate radius blows up. 13 ' 1 ^ This 
suggests that the boundary circumference is an interesting parameter to include in the 
analysis of cylinder solutions. 

In this paper we consider a circumference calculated using a Riemann-Liouville fractional 
integral. This introduces a new parameter into the cylinder description, the fractional index 
a. The circumference calculated from the Levi-Civita metric involves the angular deficit 
factor. A fractional circumference can relate the physical behavior of the angular deficit to 
a mathematical framework through b(a). Matches to an interior will provide a fractional 
parametrization of all of the cylinder parameters. An unexpected result of the fractional 
extension is a new family of solutions occuring at one of the angular deficit embedding^' 
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transition points. 

In the next section we briefly review the literature leading to limits on the range o 
and develop the fractional circumference. The idea is applied to the a = Gott-Hiscock 
constant density i-sjring solution, providing a parametrization of A. The method is extended 
to a a ^ shell^ in the third part of the paper and the behavior of A and a are linked. 
An Appendix introducing some of the basic fractional integral definitions is included. 



II. CYLINDER SOLUTIONS 



The range of a 



The Riemann tensor of metric ([I]) is zero for o = (0, 1/2) and is singular as r approaches 
zero 
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The literature contains a number of discussions that provide interior matter distributions 
motivating some of this behavior. An infinite matter cylinder is often considered because, in 
the small a limit, a test particle at rest experiences an acceleration r = —2a/r. This is t|h^ 
Newtonian acceleration for a particle a distance r from a line source of mass/length (a). 
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In addition, through metric matching, a can often be related to the matter stress-energy 



coming from the field 
Gott-Hiscock solutions^ 



nations. The identification is not absolute, r^he constant density 

n 

■A and the U(l) string solutions of Garfinkle- 18 ! both use the a = 
Levi-Civita vacuum. The exact cylindrical solid matter solutions that can be matched to 
the Levi-Civita vacuum indicate a restricted range for a, < o < 1/2. Some of the 3-dim 
matter solutions that can be matched to cr = also show restricted ranges in the cylinder 
mass/length, A, following from the field equations. Israelrshells bounding vacuum Levi-Civita 
and a second vacuum interior indicate a broader range. 

One of the early investigations of static cylindrical solutions, due to Marder^, directly 
examined the relation between a and jfehe cylinder A by matching to a Levi-Civita vacuum 
with no angular deficit. Krori and Paul^ used Marder's solution to establish a limit a < 1/2. 
For general crc-ipavidson^ developed a set of static cylindrical solutions which were used by 
Bonnor et al^^ to study the possible range of a. One of the motivating questions for 
this discussion was the zero Riemann tensors for a = (0, 1/2). The a = flat spacetime 
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was interpreted as a matter-free vacuum but the a = 1/2 case was left an open question. 
Philbin^ extended cylindrical solutions into the negative a region and suggested that the 
endpoints \a\ = 1/2 described planar rather than cylindrical matter distributions. Based on 
the behavior of gyros orbiting a cylindrical mass distribution, Herrera, Ruifernandez, and 
Santos^ also suggested that a = 1/2 described a planar mass. In a later paper Herrera 
et al treated the a = 1/2 transition point in greater detail and, with the (</?, z) coordinates 
taken as interchangeable to include possible planar topologies, suggested models for the 
coordinate range < o < oo. 

Some of the work on static cylinders is embedded in discussions of rotating cylinders 
with the parameter and density limits appearing for both cases. For example, one of the 
early indications of a restricted relativistic density range was the rotating dust cylinder of 
Vishveshwara and Winicour.l^ This interesting paper contains an expression for A related 
to an angular deficit and identifies the 6 = deficit factor as a critical limit related to a 
rotating column with A < 1/4. Building on this work, Lathrop and Orsene^ considered a 
cylinder with two counter-rotating dust currents and, for this matter source^ dupicated the 
Vishveshwara- WinicourW density limit, as did later work by da Silva et al.^' 



Fractional Circumference 



A fractional integral is a function convolution over a range. The fraet- 



to calculate the fractional circumference is the Riemann-Liouville form 
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where a is the fractional order, Re(a) > 0. A fractional circumference is calculated by 
integrating y/g^, around the circle, x is identified as <ft with the circumference following in 
the cf) — > 2tt limit. For the base metric, consider a 3+1 Minkowski metric with cylindrical 
coordinates (t,r,<fi,z). 

ds 2 = -dt 2 + dr 2 + r 2 d(p 2 + dz 2 . 

The usual circumference of a circle in the (r, 0) plane for r = tq is C = 27rro- tq is the 
coordinate radius. A fractional circumference for the same coordinate radius is 

<t> (27r) a ~ 1 

JVg^(4>-yT- 1 dy = 2nr y - 1 
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1 gives the usual circumference. For < a < 1, the circumference risi 
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usual 2irr and for a > 1 it is larger. The methods of fractional calculus^ 
very successful iiw»e|deling transport processes with anomalous microscopic time and/or 
spatial structure. 124123 One could expect that the density of matter distributions interior to 
the boundary would reflect the transport processes responsible for their growth. In this case, 
it is the variation in the circumference with a that we wish to link to matter distributions 
with an angular deficit, with the fractional variation in circumference reflecting a fractional 
surface matter distribution. If coordinate ranges are not imposed, the angular deficit can 



be transformed away. Here it is strongly he 



imposed angular coordinate ranges. Bonner 



ked to a physical description with the usual 



231 noted that the angular deficit parameter, b, 



determines the topology of the manifold covered by metric ([[]) and cannot be removed by 
scale transformations. Angular deficit is a topological defect like the gravitational Aharonov- 
Bohm effect discussed by Jensen and Kucera^ 

In the following, we apply the fractional circumference to the parameters of the a = 
Gott-Hiscock constant density string. 



The Gott-Hiscock String 

The Gott-Hiscock solution describes an interior constant density, 5, cylinder matched to 
a a = Levi-Civita vacuum with angular deficit b. The matter (-) and vacuum (+) metrics 
are 

ds l = _ c 2 dt 2 + d() 2 + [ !! n M]2 # 2 + dz 2 (4) 
ds l = _ c 2 dt 2 + dr 2 + b 2 r 2 d(j) 2 + dz 2 (5) 

The string radii in the interior and exterior are po and tq. The matching relations are 
sin(5p ) = 8br , and b = cos(<5po) and the mass per unit length, A, of the constant density 
string is A = (1— 6)/4. The angular deficit /excess associated with b is A0 = 27r(l — b) = 8ir\. 
The allowed range of 6 is— 1 < 6 < 1 and the positive A range is < A < 1/2. 

The boundary match provides an expression for the coordinate radius, r and associated 
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FIG. 1: embedding diagrams for the GH string 

circumference 

r = tan( ^ o) (6) 
C = 2nr b = 2.^1 (7) 

A (r, 0) cross section of the exterior vacuum is regarded as a circle with a missing pie slice. 
r b is the radius of the equivalent closed circle. The matching relations indicate that r 
approaches oo as Spo — > tt/2, while 6 = and the circumference takes its maximum value. 
Gott^ has motivated this behavior by embedding the (t, z) = const Levi-Civita vacuum 
into a 3-dim metric with coordinates (w, r' = br, 0), the embedding relation between w and 
r' defining a cone in the 3-dim space. This embedding picture is illustrated in the sketches of 
Figure 1. The shaded region represents a matter cap in the interior. The sketches illustrate 
the entire b range. In the integer model, as the density increases from A = (6 = 1) to 
A = 1/4 the upper sides of the cone tilt outward to form a b = cylinder. As the mass 
continues to increase, the center moves out into the vacuum, the lower sides tilting in to 
form an upward pointing cone. Conventionally this is associated with negative b, angular 
excess and a quasi-regular singularity^ in the vacuum. 

The circle defined by the matter/vacuum match has an effective radius rob and is clearly 
maximum for the 6 = cylinder. 

Fractional Angular Deficit 



The metric circumference of the string in the exterior spacetime is 2irbro where bro is the 
effective radius of the closed circle. Equating this to the a > fractional circumference one 
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FIG. 2: b vs a, a > -0.5 
finds a fractional angular deficit factor 



2 *» = 2 "'»frfr< (8) 



r(i + «)" 

a = 1 can be identified with 6=1, the vacuum case. Fig.2 shows 6(a) over the range 
— 0.5 < a < 1. 6 decreases toward as a — > and A moves toward the cone — > cylinder 
transition value, A = 1/4. 

In this region, 6r , the effective radius in the circumference, increases from its vacuum 
to its maximum value 6ro = 1/5. This behavior is implicit in the embedding diagrams. 
The value of a follows from 5po = sin _1 (6r ). In the a region shown in Fig.2, the fractional 
index is only another parameter that provides a mathematical base for 6. However, while 
the fractional circumference is strictly defined for a > 0, the graph shows a continuation 
across the origin into the undefined negative alpha region. The interpretation of negative a 
will come from the string models it parametrizes. In Fig.2, there is no identification of an a 
parametrization for 6=0 until the alpha range is extended further along the negative axis. 
This extension is shown in Fig. 3. There are a sequence of points of decreasing amplitude, 
oscillating around 6 = 0. The fractional string model has only a small extension into the 
negative 6 region with its vacuum quasi-regular singularity. 

The identification of a fractional form for 6 can be extended to the string mass/length 
showing A oscillations about the cone — > cylinder critical value A = 1/4. 
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FIG. 3: b vs a, -0.5 < a 



FIG. 4: A vs a, -0.5 < a 

The mass/length is shown in Figures 4 and 5. Fig.4 supports the identification of b = 1, 
a = 1 with a zero string vacuum and Fig. 5 shows the expected oscillations in the string 
mass/length about A = 1/4. The string coordinate radius, tq = tan(5po)/S, blows up at 
5po = tt/2 but the product br entering into the circumference is finite, showing small 
oscillations around r b = 1/5, the effective radius for the cylindrical embedding. This 
example used a 3-dim matter distribution inside the cylinder. The next example considers 
a cylindrical shell. 
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FIG. 5: Xvs a, a < -0.75 
III. CYLINDRICAL SHELLS 

2+1 cylindrical shells have also been used to explore the static parameter /dsssity limits. 
A sequence of cylindrical shells was first suggested by Marder^ and Wang et al^" discussed 
a shell with stress-energy resulting from a jump from exterior Levi-Civita to an interior 
Minkowski. A simple cylindrical model with non zero a has exterior (+) and interior (-) 
metrics 

ds\ = -(r/r ) 4a dt 2 + (r/r )-^[(r/r ) Sa2 (dr 2 + dz 2 ) + b 2 r 2 d(j) 2 } (9) 
ds 2 _ = -dt 2 + dp 2 + dz 2 + p 2 d(p 2 (10) 

The metric match across a surface (r ,p ) is 

Po = br (11) 

This matching relation suggests that the circumference is bounded as b and r change. The 
extrinsic curvatures are K^ b = n^. b = —T ab nf = d r g a b/2 at r = r . The shell has stress- 
energy, S a b, related to the extrinsic curvature jump < K ab > = K^ b — K~ b , < K > = < >, 



the boundary. ^ — ' Calculating the stress-energy one finds, with F(cr) = 1 — 4cr + 4o" 2 



across 



SirS* = < K) > -h) <K> (12) 

-8nSt = 1 ~ b / {(T) (13) 
br 

-SirS* = - 4a 2 /r (14) 

-8nS z z = (15) 

Using the same simple density integral to calculate A as was used in the Gott-Hiscock string, 
one finds the stress-energy parameters 

4A = 1-6F(<t) (16) 

8vrP = 4a 2 /r (17) 

8nP z = (18) 
br 



i' 



1 - bF(a) . . 

8tie = -— ^ 19 

br n 



The general equation of stated ^ is 

e + P^ + P z = Aa/r (20) 

Looking first at a = 0, the usual stiff string equation of state is found. This is an interesting 
comparison case to the GH string. It does describe a constant density object with an axial 
tension and A but the limits are set by the angular deficit factor, not by trignometric limits. 
Using ad hoc values b = (1,0, —1), we have 4A = (0, 1,2). Thus we have the same set of 
boundary values as the GH string. With b = (2ir) a ~ 1 /T(l + a) the mass/length has the same 
form as for the GH string 

4A = 1-^; (21) 

r(i + a) v ; 

and the same density oscillations around 6 = are present. A large difference between the 
two cases is that b is not limited to a unit range, without the GH identification of b. The 
current limits on a would restrict b. For a^0 the immediate observations are, that for fixed 
radius, a determines the tangential stress while the axial stress is determined by b. Links 
between a and b follow from assumptions on the partial equations of state linking e and 
individual stresses, or on relations between the stresses. For example, a shell with uniform 
stress, = P z imposes b = 1/(1 — 4a 2 ). In this model, isotropic stress is associated 
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with negative mass/length. A physical solution will have anisotropic stresses, avoiding a 
singularity at a = 1/2. 

IV. CONCLUSION 

In this paper we have examined a fractional parametrization b(a) of the angular deficit 
factor. The original motivation was to match the physical behavior of the angular deficit to 
a mathematical framework providing a fractional parametrization of cylinder parameters. 
For the Gott-Hiscock cylinder, we found new behavior when the angular deficit factor is 
continued from a vacuum with a conical defect across the embedding cylinder, 6 = 0, into 



y have implications for the 
This is under investigation. 



a vacuum with a quasi-regular singularity. The behavio 
stability of the cylinder solutions in the b parameter range. 
The models considered in the paper are for the Levi-Civita a in the range < a < 1/2. 
The applicability of a fractional angular deficit as developed, here, depends on an angular 
coordinate with an associated circumference. Herrera et al^ used coordinates (cp, z) with 
ranges — oo to oo, imposing a periodicity, ip = ip + 2ir on the angular coordinate for a < 1/2. 
Our result applies in that region. Their results seem to imply that a = 1/2 is associated 
with a planar topology and that there is no appropriate angular coordinate. For a > 1/2 
they suggest that the Levi-Civita coordinates (</?, z) are interchanged, with z becoming the 
periodic measure, z = z + zq. The circumference for this a region would be 



r(a)^°Q r(l + a) 

It is also possible that both z and (p are periodic in this region suggesting a toroidal matter 
distribution rather than a cylindrical one. This extended view would allow two fractional 
circumferences to be included in the modeling. The complete topological description of 
Levi-Civita matched matter distributions, and its fractional extensions, still have many 
open questions. The infinite cylinder, both static and rotating, continues to be an inter- 
esting and useful structure in general-relativity: The parameter effects due to including a 

1 J 

ginning to be examined^ and, with the 



cosmological constant in Levi-Civita^'^ ar 
recent connections to braneworld cosmologies 
physical relevance. 




relativistic cylinders have acquired new 
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the development o£-» 
Riemann and Weyl^' 



Appendix A: FRACTIONAL CALCULUS 

Fractional calculus was invented in 1695 when L'Hopital^'3 asked Leibniz about the 
meaning of his notation -r-Jf for n = 1/2. This question, also asked by Bernoulli defined the 
idea of a fractional derivative and L'Hopital's question about derivatives has developed into 
a fractional calculus, a framework for integrals and derivatives of non-integer order. Euler 
made the first real contribution to the development of fractional calculus with his 1738 
proof that the fractional derivative of x n was meaningful.^" The-4irst contribution using 
the idea is credited to Abel for his 1823 work on the tautochrone.^" The idea of fractional 
calculus has attracted the attention of many of the same people that are associated with 

,ysics and general relativity: Lagrange, Laplace, Fourier, Liouville, 
* among others. There are several definitions of fractional integrals 
depending on the integral limits. The Riemann form integrates from c to x, The Liouville 
version from — oo to x, and in this paper the Riemann- Liouville form is used from to x.^' 

r a [f(x)} = ^-Jf{y){x - yT-'dy (Al) 
1 (a)o 

a is the fractional order, Re(a) > 0. There is also a Weyl form for the fractional integral 
with a +oo limit. 

In order for the Riemann-Liouville integral to converge it is necessary that f{x) satisfy 
the condition 

f( x - 1 ) = 0{x 1 - 5 ) 5 > (A2) 

Functions obeying this condition are Riemann class; for example, the function x b , b > — 1 
is a Riemann class function. The Weyl integral is useful for functions with the convergence 
property 

f(- x ) = 0{x~ v - 5 ) 5 > 0, x -> oo (A3) 

Functions satisfying this condition are said to be of Liouville class; for example, x~ b with 
b > v > 0, is a Liouville function. The allowable parameter range eliminates constants from 
this class. There can be overlap between classes. 

The evaluation of fractional integrals uses the beta and gamma functions. The gamma 
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function is denned by 

oo 

r W := Je-n x > (A4) 



r(x + l)=xF(x), x>0 

For positive integer n, T(n + 1) = n\. T(x) is extended to negative x by the functional 
equations 

T(x) := T ( x + 1 \ -1 <X< (A5) 

Ob 

T{x) ■■= T( ; X + 2 :l -2<x<-l 
v ' x(x + 1) 

_|_ ^ j 

r(x) := — -t-^ — , -n<x<-n + l 

v 7 x(a; + l) • • • (x + n- 1) 

The beta function is defined as 

i 



B(z,w) = J x z ~\l-x) w - L dx (A6) 
o 

It is related to the gamma function 

„. . T(z)T(w) , . „. 

B(2 - w) =rbT^ (A7) 

For example, the fractional integral of /(x) = 1 is 

a; 

Do v [i] = ^f(x-y) (y - 1) dy (A8) 



Comparing with the definitions of the beta and gamma functions, we see the fractional 
integral of 1 is given by 

= r^TIf" (A9) 

Another example is the fractional integral of f(x) = x a . Using the definitions the integral 
becomes 

™ = f^TTi)^ < A10 > 

For v = 1, this becomes x a+1 /(a + 1), the usual integral of x a . For a = —1, the gamma 
function in the numerator is singular so the fractional integral of \jx is not defined. The 
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analytic continuation used in the paper assigned meaning to the singular points of the gamma 
function. 
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